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Abstract —Conventional electric field integral equation based 
theory is susceptible to the spurious internal resonance problem 
when the characteristic modes of closed perfectly conducting 
objects are computed iteratively. In this paper, we present a 
combined field integral equation based theory to remove the 
difficulty of internal resonances in characteristic mode analysis. 
The electric and magnetic field integral operators are shown 
to share a common set of non-trivial characteristic pairs (values 
and modes), leading to a generalized eigenvalue problem which is 
immune to the internal resonance corruption. Numerical results 
are presented to validate the proposed formulation. This work 
may offer efficient solutions to characteristic mode analysis which 
involves electrically large closed surfaces. 

Index Terms —Characteristic mode, combined field integral 
equation, closed surface. 

I. Introduction 

HARACTERISTIC modes (CM) are originally defined 
by Garbacz as a basis to diagonalize the scattering and 
perturbation matrices of a conducting object m The theory 
was later refined by Harrington and Mautz 12 0, where an 
electric field integral equation (EFIE) based derivation was 
formulated. The EFIE based theory for characteristic mode 
analysis becomes most widely accepted due to its math¬ 
ematical elegance, reliability and convenience for practical 
applications. In a, Nalbantoglu adopted the magnetic field 
integral equation (MFIE) to determine the characteristic modes 
of simple bodies of revolution. Unfortunately, this work has 
not drawn much attention as it appears to provide no advantage 
over the EFIE based approach. 

Theory of Characteristic mode (TCM) was recently pop¬ 
ularized in the antenna community by Cabedo-Fabres m as 
it has useful applications in antenna shape synthesis, input 
impedance matching, and radar cross-section control, etc. Al¬ 
though TCM has been shown promising for systematic antenna 
design, the theoretical and numerical aspects are relatively less 
addressed. Characteristic mode analysis requires one to solve 
a generalized eigenvalue problem as XJ n = A n RJ n , where 
the dense matrix pair (X,R) is generated from the EFIE 
impedance matrix Z. When the problem scale is small, one can 
use Schur decomposition to find all eigenpairs (eigenvalues 
and eigenvectors) of the generalized eigenvalue problem with 
a computational complexity of 0(N 3 ), where N is the number 
of unknowns. This becomes increasingly unaffordable as the 
problem scale increases. However, in many applications, only 
J n with small |A n | are desired. Iterative eigensolvers such 
as the Lanczos and Amoldi methods can be used for this 
purpose. For instance, the commercial software FEKO incor¬ 
porates the ARnoldi PACKage (ARPACK) which implements 


the Implicitly Restarted Arnoldi Method (IRAM) to compute 
a few characteristic modes of interest of more complicated 
geometries ( 30 - Such iterative methods are normally more 
efficient than the direct eigen-decomposition, since they only 
require a number of matrix-vector multiplication (MVM) 
operations, each of which has a computational complexity 
of 0(N 2 ). The complexity of each MVM can be further 
reduced to 0(N log TV) if a modified multilevel fast multipole 
algorithm (MLFMA) is employed [9jr 11 ]. 

Usually, one can easily obtain good approximations to the 
eigenvalues of largest magnitude with Krylov methods such 
as IRAM (7). Fast convergence to the desirable spectrum 
requires one to transform the aforementioned generalized 
eigenvalue problem to a standard one as X -1 RJ n = A“ 1 J n . 
Hence, a number of X -1 u has to be computed, where u are 
arbitrary vectors. In many cases, X -1 u can only be computed 
iteratively since the direct decomposition of X is not feasible. 
An equivalent eigenvalue problem for computing characteristic 
modes is given by Z -1 RJ n = (l + iA n ) -1 J n , which calls for 
iterative solutions of EFIEs Z -1 u. Apparently, it is favored to 
solve Z _1 u as there exists efficient techniques to precondition 
EFIEs. However, if the EFIE based TCM is applied to closed 
conducting surfaces, Z becomes ill-conditioned when the 
operating frequency is selected close to frequencies of spurious 
internal resonances. The internal resonance problem is even 
more severe for electrically-large objects. 

The difficulty of internal resonances has been well addressed 
in many radiation and scattering problems, where the com¬ 
bined field integral equation (CFIE) is an effective remedy 
fl2l . Although a full-rank Z can be obtained by combing the 
electric and magnetic field integral operators, it is not trans¬ 
parent to update the vector u consistently in the eigenanalysis 
where characteristic modes are iteratively computed. In this 
paper, an MFIE based TCM is formulated, which is inspired 
by the very brief result documented in (H. Then, a CFIE based 
TCM can be formulated as the EFIE and MFIE based TCMs 
share the same set of non-trivial characteristic pairs (trivial and 
spurious ones form the null space of the impedance matrix). 
Moreover, the corresponding characteristic mode expansion of 
any excited surface current is obtained. Since EFIE can be 
preconditioned with a Calderon multiplicative preconditioner 
(CMP) Mm, we can further formulate a CMP-CFIE based 
TCM, which is free from internal resonances, and can be easily 
solved using iterative eigensolvers after being transformed to 
a standard eigenvalue problem. 
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II. EFIE Based TCM 

The EFIE for an arbitrarily shaped perfect electric conductor 
(PEC) object relates the tangential incident field E inc and the 
scattered field E sca on the PEC surface S as 

nxE sco (r) = 2 E (r,r')-J(r') = -nxE itlc (r), r G S (1) 

where 

Z E (r, r') • J(r') = h x C E { r, r') • J(r') (2) 

In the above, we denote 

£ E (j, r ') • J(r 7 ) *= ikrj f dr' G(r, r') • J(r'), (3) 

Js 

and J(r') is the induced surface electric current, n is the unit 
normal of S. The Green’s dyadic is 

G(r,r') = 

with the unit dyad denoted as I. The scalar Green’s function 
is 

p ikR 

5( r , r ') = ^;, i?=|r —r'| (5) 

where r and r' denote the field and source points, respec¬ 
tively. Moreover, the wavenumber is k = while the 

wave impedance is 77 = yj y i/e , where the permittivity and 
permeability of the homogeneous medium are denoted by e 
and /I, respectively. 

The most widely adopted TCM is formulated on top of 
EFIE. As suggested by Harrington and Mautz 0, the char¬ 
acteristic modes of arbitrary PEC objects can be found by 
solving the operator eigenvalue problem given by 

X E (r, r') • J„(r') = A„7£ s (r, r') • J n (r'), r e S (6) 

where 

7 Z e = 1 -(Z e + Z* e ), X e = j.{Z e - Z* e ) (7) 

with the complex conjugate denoted by *, and J n (r') denotes 
the characteristic currents, and A n denotes the characteristic 
values. An equivalent TCM to is easily obtained as 

Ze(t, r') ■ Jn(r') = (1 + iA n )7l E (r, r') ■ J„(r'), r G S (8) 


VV 

k 2 


s(r, r') 


( 4 ) 


Characteristic values A n is important as |A n | indicates the 
modal behavior. When A n = 0, the corresponding J n is an 
externally resonant mode which is efficient in radiating energy. 
When A n > 0 (A n < 0), J n is an inductive (capacitive) mode 
which stores predominantly magnetic (electric) energy. When 
|A n | = oc, J n corresponds to the internal resonance (trivial 
and spurious) modes of the closed PEC surface {Ze • J n = 0), 
which has no contribution to radiated or scattered fields. In 
this case, Z#, R# and X# are rank-deficient, and it is not 
feasible to convert generalized eigenvalue problems and 
© to standard ones as 


X £; 1 R£;J n = (A n ) 1 J n 

(11) 

Z £ ; 1 R£;J n = (1 + i\ n ) 1 Jn 

(12) 


whose desired spectra can be easily found iteratively. 

TCM may be formulated using Yaghjian’s augmented EFIE 
(AEFIE) lITbl . By demanding that the normal component of 
the electric flux density D(r) = e:E(r) be equal to the surface 
charge density on S, or 

n ■ E(r) — t — V • J(r) = 0, r G S' (13) 

IUJ£ 

one has 

£ E (r,r') • J(r') + • J(r) = -E inc (r), reS (14) 

Taking nx and n- on both sides of ( [14] ) yield two equations 
with consistent solutions. The resultant matrix equations form 
an over-determined system whose solution can be solved using 
the method of least squares. This scheme is immune to the 
internal resonance corruption. Such AEFIE based TCM may 
be constructed as 

3to[£ e ( r, r')] • J„(r') + • J„(r) 

k ( 15 ) 

= A ra 5fte[£ B (r,r')] • J„(r') 

with r G S, which can be solved using a similar discretization 
method. Also of note is that this scheme may fail when a 
sphere is studied. 


III. MFIE Based TCM 


Approximating the characteristic currents with expansion 
functions f ? -(r / ), and weighting the two operator eigenvalue 
problems ^ and ^ with testing functions G(r), one can 
obtain two matrix eigenvalue equations as 

X#J n = A n R^J n (9) 

and 

Z^J n = (l + iA n )R^J n (10) 

In the above, R# and X# are taken as the real and imaginary 
parts of the EFIE impedance matrix Ze, respectively, while 
Ze is the matrix representation of the operator Ze- The 
eigenvectors J n contain basis expansion coefficients as entries. 
For EFIEs, Rao-Wilton-Glisson (RWG) and nxRWG basis 
functions are chosen as f ) and G, respectively. 


To obtain an MFIE based TCM, we first consider the 
characteristic field E n radiated by a characteristic current J n , 
that is 

V x V x E n — k 2 E n = ikr]J n (16) 

According to Harrington and Mautz 0 , the scattering operator 
which operates on incoming waves to yield outgoing waves 
is diagonalized when characteristic fields E n are chosen as 
the basis of outgoing waves, and their complex conjugates 
E * as the basis of incoming waves. Based on that J n is 
real (equiphase), one can easily argue that E n + E* is a 
source-free field using ( p~6]> . Hence, if source-free fields are 
written in a linear superposition form, i.e., a n (E n + E*), 
due to linearity, it is sufficient to consider a fictitious single 
mode incident field E i nc = E n + E* which impinges upon 
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a closed PEC surface. The induced surface current J s due to 
this incident field is determined by 

Z E (r,r') ■ Js(r') = -n x [E„(r) + E*(r)] 

= -nx [C E { r, r') + C* E { r, r')] • J n (r') 
=-27£ B (r,r') • J n (r'), r eS 

(17) 

Due to Equation <[8|», it is obvious that 

J s = ( 18 ) 

1 +2A n 

One can also investigate the same problem with MFIE. It 
is easy to argue that the outgoing E n corresponds to H n and 
the incoming E* corresponds to — H* by using 

H = x E (19) 

kr\ 

Therefore, the corresponding incident magnetic field is given 
by H i n c = H n — H* , which is also source-free. The induced 
surface current J s is determined by 

Z H (r,v')-3 s (r')=l t -U inc (r) 

= It • [H„(r) — H* (r)], r e5 + 

where S+ represents a surface which is infinitesimally larger 
than S, 1 1 = —n xnx = I — fin extracts the tangential field 
components, and 


Zh{ r, r) • J(r') = J(r)xft-I r JC H ( r, r') • J(r') (21) 

with 


/Cij-(r, r') • J(r') = V x f dr'G(r, r') • J(r') 

Js 

= [ dr'\7g(r,r') x J(r') 

Js 


When r G S, one has 
Zh{ r,r') • J(r') = x h 


( 22 ) 


- It P- 


‘.V. [ dr'Vg(r,r') x J(r') 

Js 


(23) 


Note that H n (r) = /C^(r, r') • J n (r'), Equation ( |2Q| ) gives 
rise to 


Z H { r, r') • J s (r') = -2 iX H (r, r') • J n (r') 


(24) 


where 


X H = Yi - Z H) 


(25) 


Since ( [20] ) and (IT] ) govern the same physical problem, the 
same surface currents are induced. Substituting ( [18] ) into ( [24] ), 
the MFIE based TCM is obtained as 


Ztf(r,r')-J n (r') = z(l+zA n )A^(r, r') • J n (r'), r G 5 (26) 
An equivalent form is given by 

r') • Jn(r') = -A n ^(r, r') • J n (r'), r G 5 (27) 

where 


Kh = \ (Zh + Z* H ) 


Normally, EFIE and MFIE share a common set of character¬ 
istic pairs (currents and values) except the null space modes. 
Note that internal resonance modes of EFIE hardly radiate, 
while those of MFIE are induced currents which radiate 
efficiently. Hence, they do not enter the common set H~2l . 

To obtain discretized eigenvalue problems, RWG and 
nxRWG basis functions are chosen as f ) and G, respectively. 
To achieve better accuracy, one can use Buffa-Christiansen 
(BC) or Chen-Wilton (CW) basis functions as expansion 
functions t*. Hence, ( [27] ) and ( [26] ) lead to 

RtfJn = — A n X^J n (29) 

and 

Zh-J n = ^(1 + ^A n )X#J n (30) 

where R h and X# are the real and imaginary parts of 
Z h, respectively, while Zh is the matrix representation of 
operator Z H . Both R h and Z h become ill-conditioned near 
frequencies of internal resonances. 


IV. CFIE Based TCM 

Following the construction of CFIE for radiation or scatter¬ 
ing problems lfl7l . the CFIE based TCM is formulated as 

[aZ E { r, r') + (1 - a)rjZ H (r, r')] • J n (r') 

= (1 + iK) [otH E { r, r') + (1 - a)irjX H (r , r')] • J n (r') 

(31) 

where r G S', and a is the combination coefficient. After 
discretization, the matrix eigenvalue equation is written as 

Z c J n = (l + iA n )K c J n (32) 

where 


Z c = o/ZiE + (1 — ol)tjZh (33a) 

K c = ckR e H - (1 — (33b) 


Since Zc is full-rank even at frequencies of internal reso¬ 
nances, (32) is transformed to a standard eigenvalue equation 
as 


Z c K(^J n — (1 + iX n ) J 


(34) 


which can be iteratively solved. 

Characteristic mode expansion of excited currents is for¬ 
mulated as follows. Consider an excitation problem which is 
determined by the CFIE as 

[aZ E (r,r’) + (1 - a)r/Z H (r,r')} -J(r') 

= -an x Ej„ c (r) + (1 - a)rjit ■ H inc (r) 
where r £ S. The discretization of ( [35] ) gives rise to 

ZcT = F inc (36) 


where I is a column vector containing the expansion coeffi¬ 
cients of RWGs f ‘j, and F inc is the vector representation of 
the mixed fields —anx E i nc + (1 — • H i nc . 

Normally, Z^ is not a symmetric matrix. One can construct 
an auxiliary eigenvalue problem of ( [32] ) as 

z£j“ = (l + U„)K£j“ 


(28) 


(37) 
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where the superscript T denotes the transpose of a matrix. The 
Kc-orthogonality is therefore obtained as 

(J°) T K c j = I (38) 

where matrices J a and J contain characteristic currents 
and J n as columns, respectively, and I is the identity matrix. 
We assume I to be a linear superposition of J n , that is 

I = Ja (39) 

where a is a vector containing the modal coefficients as entries. 
Using the orthogonality condition, then 

(J°) T Z c j = s (40) 

where the diagonal matrix X = diag[l -MAi, I-HA 2 , • • • ], one 
can obtain the modal coefficients as 

a= (S) _1 (J a ) T F inc (41) 

Taking a = 1, the above is reduced to the conventional EFIE 
case where Zc becomes symmetric. 


V. CMP-CFIE Based TCM 


In scattering or radiation problems, a Calderon multiplica¬ 
tive preconditioner (CMP) has been used to accelerate the 
convergence of CFIE solutions lfl4l fl5l . Similarly, we can 
obtain the CMP-CFIE based TCM as 


a 


-Z E (ik) • Z E (k) + (1 - a)rjZ H (k) 

T] 


— (1 +^A n ) 


a 


—ZE(ik) • 7 Z E (k) + (1 — 

T] 


(42) 

The discretization of the above simply follows lfl4lfl5l , which 
leads to 

ZccJn = (1 + i\ n )K CC Jn (43) 


where 

Zee = —ZeG~ 1 Ze + (1 — a)rjZiH (44a) 

T] 

Kcc = -ZbG-'Rb + (1 - a)ir,X H (44b) 
V 

In the above, Ze is the matrix representation of operator 
ZE{ik) with BCs and nxBCs chosen as the expansion and 
testing functions, respectively, and G is the Gramian matrix 
linking nxRWGs and BCs tl3ll . It has been shown that Zee 
is well-conditioned even at frequencies of internal resonances. 


VI. Numerical Results 

In this study, the discretized standard eigenvalue prob¬ 
lems are solved with the popular IRAM embedded in MAT- 
FAB R2014b on an Intel Core i7-4700MQ CPU with 2.40 
GHz clock rate. In each IRAM iteration, we compute the 
MVM Z -1 u with the generalized minimal residual (GMRES) 
method, where Z can be Z E , Z#, Zc. Zee , etc. In each 
GMRES iteration, an MVM in the form of Zu is directly 
computed for convenience, where Z is explicitly stored. For 
large-scale applications, one can apply MFFMA to compute 
MVMs such as Zu, Ru and Xu implicitly with a complexity 



Fig. 1. Condition numbers of EFIE, MFIE and CFIE matrices. It clearly 
indicates that CFIE removes the corruption due to internal resonance modes. 


of 0(N log N) CD. In the following examples, the total 
number of IRAM (outer) iterations is denoted as N out , and 
the average number of GMRES (inner) iterations in each 
IRAM (outer) iteration is denoted as N in . The GMRES error 
tolerance is set to 10“ 10 unless specified otherwise. 

We first consider a perfectly conducting sphere with a radius 
of 1 m where 2 280 unknowns are used. The condition numbers 
of matrices Z E , X#, Z#, R E and Zc are plotted in Fig. [l] 
where the frequency is swept at a step of 1 MHz. The spikes in 
the condition numbers of Z e correspond to spurious internal 
resonances of EFIE. It is obvious that Ze and X# have the 
same condition numbers at most frequencies except that X# 
seems to have multiple spikes around each internal resonance 
frequency. This renders it difficult to compute characteristic 
modes using © and ©• Although MFIE is usually well- 
conditioned, it also suffers from internal resonances as several 
spikes are observed in the condition numbers of Z#. So does 
R h except that it is more poorly-conditioned. Hence, (29) and 
(30) may not be suitable for characteristic mode computation 
when the operating frequency is near internal resonances. 
On the other hand, CFIE is free from the internal resonance 
corruption as no spikes are observed in the condition numbers 
of Ze within the entire frequency band. By setting the CFIE 
combination coefficient a to 0.5, Ze has a small condition 
number ranging from 8 to 17. 

A finer frequency sweep is performed with respect to the 
spikes of Z E . Two lowest frequencies of internal resonances 
are located with better accuracies, which are listed in Table |T| 
with the condition numbers of Z E , X#, Z#, R h and Ze 
provided. It clearly shows that Ze is immune to the internal 
resonance corruption which the other matrices are susceptible 
to. At the internal resonance frequency 131.192 MHz, we 
excite the sphere by a ^-polarized plane wave propagating 
along ^-direction. The induced current computed by EFIE, 
as shown in Fig. |2(a)| is incorrect as it is dominated by 
a spurious internal resonance mode. By solving a standard 
eigenvalue problem Z#u n = £ n u n , one can easily find the 
internal resonance mode which corresponds to the smallest 
|Cn| (around zero). The induced current can be correctly solved 
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TABLE I 

Condition Numbers of EFIE, MFIE, and CFIE Matrices at 
Internal Resonances (Sphere) 


Freq (MHz) 

Z E 

x £ 

Z H 

R H 

z c 

131.192 

1.02 x 10 6 

1.06 x 10® 

839 

2 437 

13 

185.005 

6.13 x 10 5 

6.44 x 10 5 

275 

388 
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for by CFIE with a combination coefficient a = 0.5, as shown 


in Fig. 2(b) 


We then compute 100 modes with the smallest |A n | by 
EFIE, MFIE and CIFE, and list the first few characteristic 
values A n in Table [H] It is observed that the discrepancies 
between the results computed by EFIE and MFIE are larger 
for the first 3 characteristic values, which correspond to MFIE 
internal resonance modes. Since MFIE is not implemented 
to generate so accurate solutions as EFIE, the numerically 
computed internal resonance frequencies of the two are not 
exactly the same, and CFIE leads to complex A n with small 
imaginary parts. Even so, the real parts of A n obtained by 
CFIE agree well with EFIE results. 

At frequencies of spurious internal resonances, EFIE may 
not be able to find correct modal currents. For example, we 
plot in Fig. 3(a) the third modal current J 3 computed by EFIE 
at 131.192 MHz. It is not consistent with those found by EFIE 
at a lower frequency, e.g., 128 MHz [Fig. |3(c)| , or a higher 
frequency, e.g., 136 MHz [Fig. 3(d)| . However, the correct 
modal current can be obtained with CFIE even at internal 
resonances, as shown in Fig. |3(b)| Furthermore, in Fig. [4| 
we plot the normalized quantities |J^Z#J n | and |J^Z#J n | 
for the first 100 modes, where J n are computed by CFIE at 
131.192 MHz. Obviously, the first 3 modes correspond to null 
space modes of Z h which are well radiating. Moreover, J 68 
(^68 = —739), J 69 (^69 — —864) and J 70 (A 70 = —2731) are 
null space modes of Z e which are poorly radiating. They have 
similar current patterns as shown in Fig. | 2 (a)| Table [Tll| shows 
AI ou t, N in and CPU times of different TCMs, where EFIE and 
EFIE* denote the cases of using (12) and (IT), respectively. 
It is obvious that CFIE based TCM is the most efficient, even 
though EFIE* involves real number evaluation only. Besides, 
in Table |!v| we compare a few characteristic values computed 
by EFIE, MFIE and CFIE when the operating frequency is 
shifted to 128 MHz. In this case, better consistency is observed 
between the results of EFIE (CFIE) and MFIE based TCMs. 

We next consider a cuboid with a dimension of 2.0 x 1.6 x 
1.2 m where 1311 edges are used. The lowest frequency of 
spurious internal resonances is found to be 119.88 MHz using 
a fine frequency sweep. By setting the operating frequency 
to 119.880 MHz, we excite the cuboid with a ^-polarized 
plane wave propagating along ^-direction. The induced current 
computed by EFIE is incorrect as it is dominated by a spurious 
internal resonance mode [Fig. |5(a)| , which again corresponds 
to a null space vector of Z e- The correct induced current 
is computed by CMP-CFIE with a combination coefficient 
a = 0.9 as demonstrated in Fig. |5(b) Table [V] lists the 
condition numbers of Z e, Z#, Zc, and Zee at an internal 



Fig. 2. Surface currents of a sphere excited by a plane wave at 131.192 MHz: 
(a) Incorrect current computed by EFIE which is dominated by a spurious 
internal resonance mode, (b) Correct induced current computed by CFIE (0.5). 

TABLE II 

Characteristic Values (CVs) Computed by EFIE, MFIE, and 
CFIE at Internal Resonance 131.192 MHz (Sphere) 


CVs 

EFIE 

MFIE 

CFIE 

Ai 

0.0503 

0.2457 

0.0503 - 0.0005* 

A2 

0.0508 

0.2697 

0.0508 - 0.0004* 

A3 

0.0509 

0.5293 

0.0509 - 0.0004* 

a 4 

1.2628 

1.2716 

1.2680 - 0.0026* 

A5 

1.2660 

1.2720 

1.2680 - 0.0031* 

A6 

1.2666 

1.2722 

1.2682 - 0.0026* 

A7 

1.2670 

1.2727 

1.2683 - 0.0022* 

As 

1.2675 

1.2734 

1.2693 - 0.0027* 

a 9 

-1.3233 

-1.3293 

-1.3277- 0.0026* 

A10 

-1.3237 

-1.3299 

-1.3280 - 0.0029* 


resonance frequency, e.g., 119.88 MHz, and a non-internal 
resonance frequency, e.g., 121 MHz, respectively. It is obvious 
that the CFIE matrix Zc is well-conditioned at the internal 
resonance frequency, while the CMP-CFIE one Zee is even 
better conditioned. This is reflected by the eigen-spectral 
distributions of Z#, Zc and Zees as plotted in Fig. [6 
At the internal resonance frequency 119.88 MHz, we com¬ 
pute 50 modes with the smallest characteristic values |A n | 
by EFIE, CFIE, and CMP-CFIE, respectively, and list the 
first few characteristic values in Table VI Good agreement is 
observed except for that A n resulted from CFIE or CMP-CFIE 
are not purely real. This is due to the accuracy inconsistency 
between EFIE and MFIE solutions. Some characteristic cur- 


TABLE III 

Iteration Numbers and CPU Times of EFIE and CFIE at Internal 
Resonance 131.192 MHz (Sphere) 


TCM 

EFIE 

EFIE* 

CFIE 

A" out 

100 

101 

100 

N in 

381.4 

391.2 

78.6 

CPU time (secs) 

789.9 

300.5 

99.9 
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Fig. 3. Sphere modal current J 3 : (a) Pattern inconsistent with others, 
computed by EFIE at 131.192 MHz. (b) Correct pattern computed by CFIE 
(0.5) at 131.192 MHz. (c) Correct pattern computed by EFIE at 128 MHz. 
(d) Correct pattern computed by EFIE at 136 MHz. The corresponding 
characteristic values of the 4 currents are 0.0509, 0.0509 — 0.0004i, 0.1100, 
and —0.0382, respectively. 
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number index 


Fig. 4. Normalized quantities |jJZ#J n | and |J^Z#J n | of the first 100 
modes where J n are computed by CFIE at 131.192 MHz. 


rents cannot be correctly computed by the EFIE based TCM 
at spurious internal resonances. We demonstrate the current 


pattern of J 3 as an example, which is shown in Fig. 7(a) 


The correct current patterns can be obtained by CFIE and 
CMP-CFIE based TCMs, which are shown in Figs. |7(b)| and 
(c), respectively. When the operating frequency is shifted off 
the internal resonance one, the correct current mode can be 
found among the EFIE results after mode tracking, as shown 
in Fig. 7(d) A comparison of the normalized far field patterns 
of this current at 119.88 MHz and 121 MHz are plotted in 
Fig. 1 8(a) and (b), respectively, where good agreement between 
the results of EFIE and CFIE is observed. Table IVIII shows 


TABLE IV 

Characteristic Values (CVs) Computed by EFIE, MFIE, and 
CFIE at 128 MHz (Sphere) 


CVs 

EFIE 

MFIE 

CFIE 

Ai 

0.1094 

0.1157 

0.1095 - 0.0004i 

A 2 

0.1099 

0.1159 

0.1099 - 0.00037 

A 3 

0.1100 

0.1183 

0.1100 - 0.00037 

a 4 

-1.2755 

-1.2814 

-1.2797 - 0.00267 

A 5 

-1.2760 

-1.2819 

-1.2799 - 0.00297 

A 6 

-1.2779 

-1.2837 

-1.2818 - 0.00277 

A 7 

-1.2795 

-1.2851 

-1.2833 - 0.00267 

As 

-1.2808 

-1.2861 

-1.2844 - 0.00257 

Ag 

1.3824 

1.3881 

1.3846 - 0.00257 

A 10 

1.3828 

1.3885 

1.3846 - 0.00327 

TABLE V 

Condition Numbers of EFIE, MFIE, CFIE and CMP-CFIE 
Matrices (Cuboid) 

Freq (MHz) 

Z E 

Z H Z c Z CC 

119.88 

7.4305 x 10 5 

106 58 10 

121 


276 

65 44 10 



(a) (b) 

Fig. 5. Surface currents of a cuboid excited by a plane wave at 119.88 MHz: 
(a) Incorrect induced current computed by EFIE which is dominated by a 
spurious internal resonance mode, (b) Correct induced current computed by 
CMP-CFIE (0.9). 


N ou t , Nin and CPU times of different TCMs. It is obvious 
that both CFIE and CMP-CFIE based TCMs are more efficient 
than the EFIE based one. The speed-accuracy trade-off needs 
to be considered as setting a large suppresses the imaginary 
parts of A n but increases the CPU time. Based on the results 
presented in [15 , 18], the CMP-CFIE scheme is promising for 
large-scale applications where the conventional EFIE fails to 
function. 

In the last example, we demonstrate the characteristic mode 
expansion with a NASA almond where 1884 edges are used. 
We first solve a scattering problem where the almond is 
illuminated by a ^-polarized plane wave propagating along 
—^-direction at 26.81 MHz. The induced current computed by 
EFIE is corrupted by an internal resonance mode, as shown 
in Fig. |9(a)| However, by setting a to 0.5, one can use CFIE 
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Fig. 6 . Eigenvalues £n of matrices Z e, Z c and Z cc at internal resonance. 
No zero eigenvalues are found for CFIE and CMP-CFIE. CMP-CFIE has the 
most compact eigen-spectral distribution. 

TABLE VI 

Characteristic Values (CVs) Computed by EFIE, CFIE, and 
CMP-CFIE at Internal Resonance 119.88 MHz (Cuboid) 


CV 

EFIE 

CFIE 

CMP-CFIE 

Ai 

0.2389 

0.2391 + 0.001U 

0.2390+ 0.0013i 

a 2 

-0.4643 

-0.4639 + 0.0022z 

-0.4637+ 0.004U 

As 

0.4890 

0.4892 + 0.0007z 

0.4891 + O.OOlOi 

a 4 

0.4979 

0.4980 + 0.0002z 

0.4979 + 0.0004i 

As 

-0.6834 

-0.6826+ 0.0038z 

-0.6823 + 0.0087z 


to find the correct induced current as illustrated in Fig. |9(b) 
We find a few important characteristic modes with the 
CFIE based TCM, and use them to reconstruct the induced 


current. Thus, only a few columns of J and J a in (41) which 
correspond to small |1 + i\ n \ are kept to compute the modal 
coefficients. The first 100 modal coefficients a n are plotted in 
Fig. 10 It shows that not all modes with large |1 + 


l-i 


can be efficiently excited as their modal-source interaction 

(j a ) T F,;, 


is negligible. Figs. 9(c)|to (f) illustrate the currents 


reconstructed using 5, 15, 25 and 75 characteristic modes, 
respectively. Figure plots RCSs (c t> = 0°) of the currents 
computed with EFIE, CFIE, as well as those reconstructed 
using different numbers of characteristic modes. It is obvious 
that the reconstructed current and the associated RCS converge 
to correct results as more modes are included. In this example, 


TABLE VII 

Iteration Numbers and CPU Times of EFIE, CFIE, and CMP-CFIE 
at Internal Resonance 119.88 MHz (Cuboid) 


TCM 

EFIE 

CFIE 

CMP-CFIE 

A" out 

100 

100 

100 

N in 

265.7 

192.3 

80.2 

CPU time (secs) 

240.0 

141.2 

44.2 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.05 0.1 0.15 0.2 0.25 0.3 0.35 


(c) (d) 

Fig. 7. Cuboid modal current J 3 : (a) Pattern inconsistent with others, 
computed by EFIE at 119.88 MHz. (b) Correct pattern computed by CFIE at 
119.88 MHz. (c) Correct pattern computed by CMP-CFIE at 119.88 MHz. (d) 
Correct pattern computed by EFIE at 121 MHz. The corresponding character¬ 
istic values of the 4 modes are 0.4890, 0.4892+ 0.0007z, 0.4891 + 0.004U, 
and 0.4705, respectively. 




(a) (b) 

Fig. 8. Far held patterns of cuboid mode J3: (a) at 119.88 MHz. (b) at 
121 MHz. 


good agreement is observed when 75 modes are considered. 
Hence, the system’s order can be greatly reduced from 1 884 
to 75, which offers a model order reduction (MOR) based on 
characteristic mode expansion. 

VII. Conclusion 

A CFIE based theory of characteristic mode (TCM) is 
presented in this paper to overcome the difficulty of slow 
convergence of EFIE solutions around frequencies of spurious 
internal resonances. Since the MFIE based formulation shares 
a common set of non-trivial characteristic pairs with the EFIE 
based one, they are combined to form a generalized eigenvalue 
problem which can be easily casted into a standard one 
where the spurious internal resonance corruption is removed. 
A CMP-CFIE based TCM is further formulated to enhance the 
performance of the proposed scheme. MOR based on the char¬ 
acteristic mode expansion is also presented which may serve 






































Fig. 11. RCSs of directly computed and reconstructed currents. 



123456 1234567 

xIO' 3 xIO' 3 


(e) (f) 

Fig. 9. Directly computed and reconstructed currents on a NASA almond 
at internal resonance 26.81 MHz: (a) EFIE solution, (b) CFIE solution, (c) 
Current reconstructed using 5 modes, (d) Current reconstructed using 15 
modes, (e) Current reconstructed using 25 modes, (f) Current reconstructed 
using 75 modes. 



number index 


Fig. 10. Modal coefficients \a n \ of the first 100 characteristic modes. 

as a useful tool in system design and optimization. The validity 
and efficiency of the proposed scheme are demonstrated in a 
few numerical experiments. It may be incorporated with fast 
algorithms such as MLFMA to compute characteristic modes 
of large-scale geometries. 
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